Finite temperature QCD in the quark-composites approach 



F. Palumbo * 
INFN - Laboratori Nazionali di Frascati 
P. 0. Box 13, 1-00044 Frascati, ITALIA 
Internet: palumbof @lnf . inf n . it 

We investigate QCD at finite temperature in the quark composites approach, which is based on 
the use of quark composites with hadronic quantum numbers as fundamentai variabies. We find 
that chiral symmetry restoration and quark deconfinement are one and the same first order phase 
transition, whose criticai temperature, in a one loop approximation, is T — 2\fQ.p~ 2 m- n , where m-n 
is the pion mass, Q. = 24 the number of up and down quark components, and p a parameter of order 
1 whose precise value can be determined by the study of the pion-pion interaction. 
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The idea that a spontaneously broken symmetry can be 
restored at sufficiently high temperature jjj has found a 
natural application in QCD. Since quantitative investiga- 
tions of this problem rely so far on numerical simulations, 
it is not easy to establish whether the chiral symmetry 
restoration is a first or second order phase transition, a 
feature which according to many authors depends on the 
number of flavors Q . The situation is somewhat compli- 
cated by the possible interplay between the restoration 
of chiral symmetry and quark deconfinement, which in a 
recent paper have been found to coincide || . 

We study QCD at finite temperature in the quark com- 
posites approach (J] with two flavors. We find that chi- 
ral symmetry restoration and quark deconfinement are 
indeed one and the same first order phase transition. Al- 
though the validity of our resul is restricted to two flavors, 
the reason why the transition is first order does not seem 
to depend on the number of flavors. We find the critical 
temperature 

T = — 5- m n . (.1) 
P- 

In the above equation is the pion mass, 17 = 24 is the 
number of up and down quark components, and p is a pa- 
rameter of order 1 whose precise value can be determined 
by the study of the pion-pion interaction. 

While the conclusion about the nature of the phase 
transition is unavoidable in the present context, we can- 
not say anything about the accuracy of the value of the 
critical temperature without an estimate of its correc- 
tions. Here it could help a numerical simulation done 
using the action of the quark composites approach. In- 
deed this approach is based on the idea that if a sig- 
nificant part of the binding mechanism of hadrons is ac- 
counted for by the use of quark composites with hadronic 
quantum numbers as fundamental variables, the residual 
interaction can be sufficiently weak for a perturbative 
treatment. If this is true the present formulation should 
also make it easier the numerical evaluation of the quark 
determinant. 

The idea at the basis of the quark composites approach 
is realized by assuming the composites with the quan- 
tum numbers of the nucleons as new integration variables 
in the Berezin integral which defines the partition func- 
tion H , and replacing the composites with the quantum 
numbers of the chiral mesons by auxiliary fields || . Then 
irrelevant operators which provide the kinetic terms for 
the composites are added to the standard action, and the 
quark action is treated as a perturbation. The expan- 
sion parameter cannot obviously be the gauge coupling 
constant, but it is instead a dimcnsionlcss constant ap- 
pearing in the definition of the composites. In conclusion 
the perturbative expansion is formulated in terms of phc- 
nomenological fields, whose effective coupling constants, 
however, are given to any order by integrals, depend- 
ing on the order, over the gluon and quark fields. The 



gluon fields can be treated, according to the dynamical 
regime, perturbatively or non perturbatively wr to the 
gauge coupling. It is an intrinsic, fundamental feature 
of the approach that the quarks, as a consequence of 
the spontaneous chiral symmetry breaking are nonprop- 
agating particles which appear only in virtual states: In 
other words their confinement is built in already at the 
perturbative level, unlike the standard perturbation the- 
ory where the quarks appear as physical particles whose 
confinement is intrinsically nonperturbative. 

The framework of our approach appears most conve- 
nient for the study of the phase transitions associated to 
chiral symmetry restoration and quark deconfinement. In 
this paper we investigate QCD at finite temperature and 
zero barion density, a purpose for which it is sufficient to 
introduce only the chiral composites. We adopt a regu- 
larization on a euclidean lattice because it seems natural 
dealing with composites, and it allows us to handle also 
the nonperturbative regime of the gluons. We assume 
the modified partition function 



Z = f[dU] J[d\d\] exp[-Sy M -S c - S Q 



(■2) 



where Sym is the Yang-Mills action, Sq is the action of 
the quark fields and Sc is an irrelevant operator which 
provides the kinetic terms for the quark composites with 
the quantum numbers of the chiral mesons. A" a is the 
quark field with color, isospin and Dirac indices a, r and 
a respectively, and the gluon field is associated to the 
link variables U^. 

All the elementary fields live in an euclidean lattice of 
spacing a, whose sites are identified by fourvectors x of 
spatial components Xk = 0, ...N and time component t = 
0, ...Nt, and satisfy periodic boundary conditions, with 
the exception of the quark fields which are antipcriodic 
in time: X(x) = X(x + Nek) = —X(x + N t et), and e* 
being the unit vectors in the k and time directions. The 
temperature is T — (iV t a) _1 . 

The quark action, with the notation 



(f,g) = a 4 J2f(x)9(x), 



can be written 



S Q = (\,Q\). 
The quark wave operator is 

Qx, v = m Q S xy + -^^l^U ^{x)5 x+e ^, y . 



We adopt the standard conventions 



(.3) 



(•4) 



■5) 



Ate {-4,..., 4}, e_ M 
U^{x) = U+{x-e ll ). 



(.6) 



2 



and omit the Wilson term for a reason to be explained 
later. 

The chiral composites are the pions and the sigma 



7r/j = i k n a 2 A75T> l A, <j = k n a 2 XX. 



(•7) 



75 is assumed hermitean and the r^'s are the Pauli ma- 
trices. An arbitrary factor of the dimension (length) 2 , 
necessary to give the chiral composites the dimension of 
a scalar field, has been written for convenience in the 
form a 2 ^. The composite action must be 0(4) invari- 
ant, so that with the inclusion of a linear breaking term 
it must have the form 



Sa = 



(• 



with 



7T? + 7T 2 + 7T 2 . 



The factor V" (we remind 
that Q is the number of quark components), has been 
introduced for later convenience. We choose |J the wave 
operator 



C 



i 2 a 



(f 



(.9) 



where □ is the laplacian on the lattice. The irrelevance 
of Sc requires that the parameter p be independent on 
the lattice spacing and that k„ do not diverge in the con- 
tinuum limit. These constraints, with the dependence on 
a assumed below for the breaking parameter m, ensures 
also the irrelevance of the chiral symmetry breaking term. 

We will now sketch the derivation of a perturbative 
expansion for the partition function of QCD. We start 
by using the Stratonovich-Hubbard transformation J7J to 
replace the chiral composites by auxiliary fields 



exp[-5c] = 



dx 



exp[— S x ] A exp ( A, —DX 



(.10) 



where 

S x = -\ P A [(X, (a 4 C)- x x) + W>, (a 4 C)"V)] 
A = JJ(detD(aj)) -1 , 

X 

D = ak^ \j> 2 4> + y/cim + ip 2 ^T ■ \ ■ 
The partition function can then be written 

dx 



(.11) 



[dU] exp[-SVM] 



A / [dXdX] exp ( A, (-D + Q)X 



exp[-S x ] 



(.12) 



The minimum of S x occurs at x = 0) but <f> ^ 0, imply- 
ing a breaking of the chiral symmetry. Since £1 is a rather 
large number we can apply the saddle-point method and 
perform an expansion in inverse powers of this parame- 
ter 0. The validity of this expansion is subject || to 
the condition p ~ 1 mentioned in the introduction. S x 
contains a quadratic part plus interactions which vanish 
in the continuum limit (but have finite quantum effects). 
The masses associated to the fields <f>, x an d the mass of 
the quarks as functions of <f> are 



1 




1 



Mq = K (p <j) + V f2m) + m,Q . 



The minimum of S x wr to 



occurs at 



a p 



(i-0, 



where 



£ = — am. 
2p 



.13) 



(.14) 



(.15) 



The mass of the pion coincides with m_. To have m_ 
Td-x we must assume the breaking parameter 



1 2 

m = —a ui, 
P * 



(.16) 



The mass of the a, which is equal to m + , diverges in the 
continuum limit. 

The quarks have a mass Mq ~ y/Qpk w /a, which seems 
also to diverge in the continuum limit. Whether this is 
actually true depends on the behaviuor of k^ with the 
lattice spacing. An exemple of the possible need of a 
scaling of k v with a is provided by the expression of the 
pion-nucleon coupling constant But irrespective of 
this dependence the quarks are nonpropagating at the 
perturbative level. Since their mass is of order yO, the 
saddle point expansion results to be a hopping expan- 
sion as far as the quarks are concerned, and as a conse- 
quence, the quark correlator connecting two sites which 
are n lattice spacings apart from one another gets the 
first nonvanishing contribution at an order of the expan- 
sion not smaller than n: The quarks are never produced 
to any finite order. In this sense we can say that they 
are confined at the perturbative level. 

In connection with the nature of the phase transition, 
it is essential to observe that the above result is solely due 
to the spontaneous chiral symmetry breaking, in whose 
absence there would be no saddle point expansion, and 
does not depend on the explicit breaking term. Since 



3 



the quark confinement is due to the spontaneous chiral 
symmetry breaking, the chiral symmetry restoration co- 
incides with the quark deconfinement: They are one and 
the same phase transition. 

We also note that since in the confining phase the 
quarks have no poles whatsoever, it seems that we do 
not have to worry about the spurious ones. In any case 
we can safely assume the quark Wilson parameter r q of 
order so that we can neglect the quark Wilson 

term to leading order in our expansion. 

Introducing the fluctuation of the 0-ficld 

9 = <t>-~$, (.17) 

we arrive at the desired expansion of the partition func- 
tion of QCD in terms of the auxiliary fields. Since only 
terms with an even power of Sq can contribute, we have 



V = -^P 2 ^ 2 ~ ln (Vttam + p 2 a<^j 



const + p 2 ip 2 . 



(.23) 



V\ is the one loop contribution at zero temperature. Its 
quadratic approximation, under the usual normalization 
conditions 



d d 2 
— Vi = Vi = 

dip oip z 



(.24) 



must vanish identically. This has been checked by an ex- 
plicit evaluation of the counterms 5m and 5p necessary to 
cancel out the one loop divergences. The calculation has 
been performed for simplicity by considering the action 
S x in a continuous euclidean space, whith a cutoff 1/a. 
Then the contribution coming from the counterterms is 



Z = / [dU] exp[-S Y M] 



d-X 



exp[— S x ] A 

/ ^ E T^v ( S Q) 2r CX P U ~ DX = E Z r- (• 

J r _Q \ )• \ J r _Q 

This is an expansion in inverse powers of k^, because 
obviously Z r ~ k^ 2r - Each term Z r is a function of ft 
which can be given a series in f2 _1 . The study of the 
pion-pion interaction determines the value of p and the 
scaling of k„ with the lattice spacing [^) . Note that we do 
not need to treat the gluon field perturbatively: If and 
where this can possibly be done remains here an open 
question. 

Finally we evaluate the critical temperature in a one 
loop approximation, by the condition that the effective 
potential be a minimum for vanishing pion mass, namely 
for 



Vi = 



b 2 a 2 



aSr 



5f_ 
~pb 



18) 



1 l\ t 2 P ( , Sp 2 

2 + W) Sp + lb 
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The total zero temperature one loop potential is 
V 1 = V 1 +F(ml) + 3F(m 2 _), 

where 



(.25) 



(.26) 



F(m 2 ) = 



1 



16(2tt) : 



1 a 9 2 4 i 

—to + 4a to — to In — — - 

2 a z m z 



(.27) 



The approximations to to^., quadratic in the field ip, are 
2 V / t p 2 3p 2 



- y/Q 

- v c 



(.19) 



2p 2 



P 2 

=a (p - 



2n 



4 2 

a ip 



P 2 . 3 P 4 2 



(.28) 



For this value of <f> the mass of the a remains divergent 
in the continuum limit: the chiral symmetry restoration- 
quark deconfinement is a first order phase transition. 
It is convenient to introduce the field ip 

0(T) = + a<^, (.20) 

because the pion mass vanishes for a finite value of it 



^From the cancellation of the divergent terms we get 
1 p* 



5p" 



5m 



(2tt) 2 n 

1 1 L 



a (2tt) 2 Q 



1 

4 

V 







3 2 


2 2 




— ln 


6 2 


2 2 



, (.29) 



if : 



2p- 



n 2 

3 TO,. 



(.21) 



Then, since ip has dimension 2, we can confine ourselves 
to approximations quadratic in this field. Let us separate 
the effective potential in its various parts 



which make V\ identically zero. 

Finally Vt is the temperature dependent one loop part 

V T = V + {T)+W-{T) - V + (0) +3V~(0), (.30) 

where 

T 



V (T) 



V = V + Vi + V T . 
Vq is the classical potential 



(.22) 



2(2tt) 3 

E ln 



— ? • 2 

a sin 



2tt 
— i 

N t 



(.31) 
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has been evaluated in M for T 3> m. Since we do 
not need its full expression, but only its derivative, we 
can avoid such limitation. Now 







V ± = 



T 



2(2tt) 3 



d 3 P 



(.32) 



Taking the limit a — -> at finite T we get 



d 



T 



+ 00 



d 3 p 



2(2tt) 3 



(.33) 



We see that the a does not contribute because of its di- 
vergent mass, while the pion contribution for vanishing 
pion mass is 



^ T [V-(T)-V-(fl)]= — 



dmt WJ 2(2tt) : 



d 3 p 



In conclusion 



n 3 



(.34) 



(.35) 



^From the condition that the effective potential be a min- 
imum for vanishing pion mass, namely for ip = Tp we get 
the critical temperature 



2Vn 



(.36) 
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To have an estimate of the accuracy of this result one 
should evaluate the two loop contribution to the effec- 
tive potential taking into account the quark contribution 
in the expans ion in fc" 1 . Another possibility is to use the 
formulation (.12) in a numerical simulation. This might 
have the advantage of an easier evaluation the quark de- 
terminant and of the guidance of the above result. 

We conclude by observing that our approach appears 
promising also in the study of the phase transition at 
high baryon density [[To| because, unlike the standard 
way, where the chemical potential acts on the quarks to 
allow the integration over the latters to be performed an- 
alytically, we can couple it directly to the nucleons, with 
a significant simplification. 
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